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511 keV line from millisecond pulsars in the Galactic center 
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Abstract. Observations of a strong and extended 
positron-electron annihilation line emission in the Galac- 
tic center (GC) region by the SPI/INTEGRAL arc chal- 
lenging to the existing models of positron sources in the 
Galaxy. In this paper, we study the possibility that pul- 
sar winds from a millisecond pulsar population in the GC 
produce the 511 keV line. Our preliminary estimations 
predict that the e^ annihilation rate in the GC is around 
5 x 10 42 s _1 , which is consistent with the present obser- 
vational constraints. Therefore, the pairs from pulsar 
winds can contribute significantly to the positron sources 
in the Galactic center region. Furthermore, since the dif- 
fusion length of positrons is short in the magnetic field, we 
predict that the intensity distribution of the annihilation 
line should follow the distribution of millisecond pulsars, 
which should then correlate to the mass distribution in 
the GC. 



1. 511 keV electron-positron annihilation 
emission in the Galaxy 

Since the first detection (Johnson & Haymes 1973) and 
subsequent identification (Leventhal et al. 1978) of the 
Galactic 511 keV annihilation line, the origin of the galac- 
tic positrons has become a lively topic of scientific de- 
bate. With the launch of the INTEGRAL gamma-ray 
observatory in 2002, the SPI board on telescope pro- 
vides a strong constraints on the morphology and inten- 
sity of the 511 keV line emission from the Galactic cen- 
ter (Knodlseder et al. 2003; Jean et al. 2003). The data 
analyses show that the line emitting source is diffuse, and 
that the line flux within 5° of the GC amounts to ~ 
(9.9±4) x 10" 4 photon enr 2 s" 1 (Knodlseder et al. 2003), 
corresponding to a luminosity of ~ 10 36 erg s _1 . The high 
line luminosity requires that the positron injection rate in 
the GC should be around (3 - 6) x 10 42 e+s -1 . Recently, 
analyses of the observational data by SPI/INTEGRAL 
with deep Galactic center region exposure show that the 
spatial distribution of 511 keV line appears centered on the 
Galactic center (bulge component), with no contribution 
from a disk component (Teegarden et al. 2005; Knodlseder 



et al. 2005; Churazov et al. 2005), and the positron injec- 
tion rate is up to 10 43 e + s _1 within ~ 6°. 

The SPI observations present a challenge to the present 
models of the origin of the galactic positrons, e.g. su- 
pcrnovae. Recently, Cassc et al. (2004) suggested that 
hypcrnovac (Type Ic supernovae/gamma-ray bursts) in 
the Galactic center may be the possible positron sources. 
Moreover, annihilations of light dark matter particles into 
e ± pairs (Boehm et al. 2004) have been proposed to be the 
potential origin of the 511 keV line in the GC. Cheng et al. 
(2006) also propose that the continuous capture of stars 
by the supermassive black hole at Sgr A* could explain 
the morphology and intensity of the e ± annihilation line. 
Here, we suggest that a population of millisecond pulsars 
(MSP) can significantly contribute to the 511 keV line in 
the GC. 



2. Motivations of a millisecond pulsar population 
in the GC 

Millisecond pulsars are old pulsars which could have been 
members of binary systems and been recycled to millisec- 
ond periods, having formed from low mass X-ray bina- 
ries in which the neutron stars accreted sufficient matter 
from cither white dwarf, evolved main sequence star or 
giant donor companions. The current population of these 
rapidly rotating neutron stars may either be single (having 
evaporated its companion) or have remained in a binary 
system. In observations, generally millisecond pulsars have 
a period < 20 ms, with the dipole magnetic field < 10 10 G. 
Figure 1 shows the distribution of observed MSPs in our 
Galaxy, and they distribute in two populations: the Galac- 
tic field (1/3) and globular clusters (2/3). In the Galactic 
bulge region, there are four globular clusters, including the 
famous Terzon 5 in which 27 new millisecond pulsars were 
discovered (Ransom et al. 2005). 

Recently, deep Chandra X-ray surveys of the Galactic 
center (GC) revealed a multitude of point X-ray sources 
ranging in luminosities from ~ 10 32 — 10 35 ergs s _1 (Wang, 
Gotthclf, & Lang 2002a) over a field covering a 2 x 0.8 
square degree band and from ~ 3 x 10 30 — 2 x 10 33 ergs 
s _1 in a deeper, but smaller field of 17' x 17' (Muno et al. 
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Observed MSPs in the Galaxy 
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Fig. 1. The distribution of the observed millisecond pul- 
sars in the Milk Way. The grey contour is the electron 
density distribution from Taylor & Cordes (1993). 

2003). More than 2000 weak unidentified X-ray sources 
were discovered in the Muno's field. The origin of these 
weak unidentified sources is still in dispute. Some source 
candidates have been proposed: cataclysmic variables, X- 
ray binaries, young stars, supernova ejecta, pulsars or pul- 
sar wind nebulae. 

EGRET on board the Compton GRO has identified a 
central (< 1°) ~ 30MeV-10 GeV continuum source (2EG 
J1746-2852) with a luminosity of ~ 10 37 erg s" 1 (Mattox 
et al. 1996). Further analysis of the EGRET data obtained 
the diffuse gamma ray spectrum in the Galactic center. 
The photon spectrum can be well represented by a broken 
power law with a break energy at ~ 2 GeV (see Figure 2, 
Mayer-Hasselwandcr ct al. 1998). Recently, Tsuchiya et al. 
(2004) have detected sub-TeV gamma-ray emission from 
the GC using the CANGAROO-II Imaging Atmospheric 
Chercnkov Telescope. Observations of the GC with the 
air Cerenkov telescope HESS (Aharonian et al. 2004) also 
have shown a significant source centered on Sgr A* above 
energies of 165 GeV. Some models, e.g. gamma-rays re- 
lated to the massive black hole, inverse Compton scatter- 
ing, and mesonic decay resulting from cosmic rays, are 
difficult to produce the hard gamma-ray spectrum with 
a sharp turnover at a few GeV. However, the gamma-ray 
spectrum toward the GC is similar with the gamma-ray 
spectrum emitted by middle-aged pulsars (e.g. Vela and 
Geminga) and millisecond pulsars (Zhang & Cheng 2003; 
Wang et al. 2005). In Figure 2, we can see that the su- 
perposed spectrum of 6000 MSPs could significantly con- 
tribute to the observed GeV spectrum (Wang et al. 2005). 



So we will argue that there possibly exists a pulsar 
population in the Galactic center region. Firstly, normal 
pulsars are not likely to be a major contributor according 
to the following arguments, the birth rate of normal pul- 
sars in the Milky Way is about 1/150 yr (Arzoumanian, 
Chernoff, & Cordes 2002). As the mass in the inner 20 
pc of the Galactic center is ~ 10 8 Mq (Launhardt, Zylka, 
& Mezger 2002), the birth rate of normal pulsars in this 
region is only 10 -3 of that in the entire Milky Way, or ~ 
1/150 000 yr. We note that the rate may be increased to 
as high as ~ 1/15000 yr in this region if the star formation 
rate in the nuclear bulge was higher than in the Galactic 
field over last 10 7 - 10 s yr (see Pfahl et al. 2002). Few 
normal pulsars are likely to remain in the Galactic cen- 
ter region since only a fraction (~ 40%) of normal pulsars 
in the low velocity component of the pulsar birth veloc- 
ity distribution (Arzoumanian et al. 2002) would remain 
within the 20 pc region of the Galactic center studied by 
Muno et al. (2003) on timescalcs of ~ 10 5 yrs. Mature 
pulsars can remain active as gamma-ray pulsars up to 10 6 
yr, and have the same gamma-ray power with millisecond 
pulsars (Cheng et al. 2004) , but according to the birth rate 
of pulsars in the GC, the number of gamma-ray mature 
pulsars is not higher than 10. 

On the other hand, there may exist a population of 
old neutron stars with low space velocities which have not 
escaped the Galactic center (Belczynski & Taam 2004). 
Such neutron stars could have been members of binary 
systems and been recycled to millisecond periods, having 
formed from low mass X-ray binaries in which the neutron 
stars accreted sufficient matter from either white dwarf, 
evolved main sequence star or giant donor companions. 
The current population of these millisecond pulsars may 
either be single or have remained in a binary system. The 
binary population synthesis in the GC (Taam 2005, pri- 
vate communication) shows more than 200 MSPs are pro- 
duced through recycle scenario and stay in the Muno's 
region. Non-thermal emissions from MSP wind nebulae 
can contribute to the X-ray sources observed by Chandra 
(Cheng, Taam, Wang 2006). 

3. Millisecond pulsars as the positron sources 

It is well known that relativistic particles from pulsar 
winds interacting with the interstellar medium form the 
synchrotron wind nebulae. Chi, Cheng, & Young (1996) 
proposed that the bulk of the comic positrons could orig- 
inate from pulsar winds. In this work, we consider that 
the electron-positron pair production occurs in the pulsar 
outer-magnetospheric region. 

It has been proposed that there is a strong multipolc 
magnetic field near the stellar surface, although a global 
dipolc magnetic field gives a good description of the mag- 
netic field far from the star (Ruderman & Sutherland 1975: 
Ruderman 1991). The typical radius of curvature I of the 
local magnetic field is on the order of the crust thickness 
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and then the total pair production rate can be estimated 
as 




ml ' ' 
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where / ~ 5.5P 26 / 21 -B 12 4 ^ 7 is the fraction size of the outer 
gap, and 



N = 2.7 x 10 27 P- 2 B CI , 9 ( y) 3 s -1 



(5) 
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Fig. 2. The diffuse gamma-ray spectrum in the Galactic 
center region within 1.5° and the 511 keV line emission 
within 6° (from Wang 2006). The solid and dashed lines 
are the simulated spectra of 6000 MSPs according to the 
different period and magnetic field distributions in globu- 
lar clusters and the Galactic field respectively. 



of the star (i.e. I ~ 10 5 cm), which is much less than the 
dipolc radius of curvature of dipole field component near 
stellar surface. The relation between the local multipole 
magnetic field and dipole field can be given by (Zhang & 
Cheng 2003) 

bA, (1) 



I 



where Bd is the dipole magnetic field of a pulsar, R is 
the radius of neutron stars. For MSPs, typically Bd ~ 
10 8 - 10 9 G, B s ~ 10 11 - 10 12 G which is much lower than 
the quantum critical magnetic field B q ~ 4.4 x 10 13 G, so 
pair cascades are also efficient in the local multipole field. 

The pair production mechanism is a synchrotron pho- 
ton cascade in a strong magnetic field. Photons will be 

converted into e ± pairs in the local magnetic field when 4. Discussions and conclusion 



is the primary electron-positrons passing through the po- 
lar gap (Goldreich & Julian 1969). Taking the typical pa- 
rameters P = 3 ms, Bd = 3 x 10 8 G, the positron injection 
rate for a MSP: N e ± ~ 5 x 10 37 e+ s" 1 (Wang et al. 2006). 

Since these pairs are created close to the stellar surface 
and the field lines arc converging, only a small fraction 
may keep moving toward the star and annihilate on the 
stellar surface. Ho (1986) showed that the loss cone for 
these pairs will approach ir/2, in other words, most pairs 
will be reflected by the magnetic mirroring effect and then 
move toward the light cylinder. These particles will flow 
out with the pulsar wind and be accelerated by the low- 
frequency electro-magnetic wave. 

Then how many MSPs in the region of annihilation 
emissions? Figure 2 has shown that 6000 MSPs can con- 
tribute to gamma-rays with 1.5°, and the diffuse 511 keV 
emission have a size ~ 6°. We do not know the distri- 
bution of MSPs in the GC, so we just scale the num- 
ber of MSPs by 6000 x (6°/1.5 ) 2 ~ 10 5 , where we as- 
sume the number density of MSPs may be distributed as 
Pmsp tx r^ 1 , where r c is the scaling size of the GC. Then 
a total positron injection rate from the millisecond pul- 
sar population is ~ 5 x 10 42 e + s _1 which is consistent 
with the present observational constraints. What's more, 
our scenario of a millisecond pulsar population as possible 
positron sources in the GC has some advantages to explain 
the diffuse morphology of 511 keV line emissions without 
the problem of the strong turbulent diffusion which is re- 
quired to diffuse all these positrons to a few hundred pc. 



their energy satisfies (Ruderman & Sutherland 1975) 
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from the outer-gap have the energy 
5.7 x 10 12 P 1 / 3 eV, so generally, the en- 
ergies of primary curvature photons and secondary syn- 
chrotron photons are higher than E CI i t , a photon-electron 
cascade will start and develop until this condition fails. 
At the end of a cascade, each incoming primary electron- 
positron can produce, on average, 



N. 



E rr \ 



1.9 x l0 3 B d , 9 P 1 / 3 (y) 3 , (3) 



In the present work, we suggest that there exists three 
possible MSP populations: globular clusters; the Galactic 
field; the Galactic Center. The population of MSPs in the 
GC is still an assumption, but it seems reasonable. Impor- 
tantly, the millisecond pulsar population in the Galactic 
center could provide the major sources of positrons. In ad- 
dition, the scenario of a MSP population in the GC can 
reasonably explain the diffuse morphology of 511 keV line 
emissions. 

Since there are many possible positron sources at 
present as noted in §1. Thus, how could we distinguish 
the model of a millisecond pulsar population from other 
models? Firstly, we can estimate the typical spatial dif- 
fusion scale of positrons in the magnetic field of the GC, 
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which is given by Xdiff ~ (Vict) 1 / 2 (Wang et al. 2006), 
where tl ~ E e /eB is the Larmor radius, E e is the energy 
of positrons, B ~ 10~ 5 G is the average magnetic field in 
the GC (LaRosa et al. 2005). The average cooling time t 
of positrons in the GC is ~ 10 6 years, so the characteristic 
diffusion scale is about 1 pc. Because of the low angular 
resolution of SPI/INTEGRAL (about 2 degrees), we can 
assume that the positrons annihilate in the local region as 
their sources, i.e. the millisecond pulsars. 

Therefore, we predict that the spatial intensity distri- 
bution of the annihilation lines should follow the spatial 
distribution of MSPs if a millisecond population exists in 
the GC. We could assume the spatial distribution of MSPs 
should follow the mass distribution of the GC though we 
do not know how well they follow each other. But because 
the proper motion velocity of MSPs is relatively low, we 
could reasonably assume that the two distributions are 
quite close to each other. Then if the positron sources orig- 
inate in the MSP population, the 511 keV annihilation line 
intensity would follow the mass (e.g. stars) distribution of 
the Galactic center region. If the positrons originate from 
supernovae or hypernovae, the 511 keV line emission could 
follow the distribution of massive stars and dense molecu- 
lar clouds. For the light dark matter scenario, the annihi- 
lation emissions may correlate to the dark matter density 
profile. Discrimination of these possible correlations may 
be tested in the future high resolution observations. 
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